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Using simple physical arguments we investigate the ca- 
pabilities of a quantum computer based on cold trapped 
ions. From the limitations imposed on such a device by 
spontaneous decay, laser phase coherence, ion heating and 
other sources of error, we derive a bound between the num- 
ber of laser interactions and the number of ions that may 
be used. The largest number which may be factored using 
a variety of species of ion is determined. 

PACS numbers: 32.80.Pj, 42.50.Vk, 89.80. +h 

In a quanttim computer binary numbers can be 
represented by quantum states of two-level systems 
("qubits"), bringing a new feature to computation: 
the ability to compute with coherent superpositions of 
numbers Because a single quantum operation can 
affect a superposition of many numbers in parallel, a 
quantum computer can efficiently solve certain classes 
of problems that are currently intractable on classical 
computers, such as the determination of the prime fac- 
tors of large numbers |^ . These problems are of such 
importance that there is now considerable interest in 
the practical implementation of a quantum computer 
1,1 . There are three principal challenges which must 
be met in the design of such a device: the qubits must 
be sufficiently isolated from the environment so that 
the coherence of the quantum states can be maintained 
throughout the computation; there must be a method 
of manipulating the states of the qubits in order to ef- 
fect the logical "gate" operations; and there must be a 
method for reading out the answer with high efficiency. 

Cirac and ZoUer |^ have made the most promising 
proposal for the implementation of a quantum com- 
puter so far. A number of identical ions are stored and 
laser cooled in a linear radio- frequency quadrupole trap 
to form a quantum register. The radio-frequency trap 
potential gives strong confinement of the ions in the Y 
and Z directions transverse to the trap axis, while an 
electrostatic potential forces the ions to oscillate in an 
effective harmonic potential in the axial direction (X). 
After laser cooling the ions become localized along the 
trap axis (the Lamb-Dicke regime) with a spacing de- 
termined by their Coulomb repulsion and the confining 
axial potential. The normal mode of the ions' collective 
oscillations which has the lowest frequency is the ax- 
ial center of mass (CM) mode, in which all the trapped 
ions oscillate together. A qubit is the electronic ground 
state |.g)(|0)) and a long-lived excited state |e)(|l)) of 



the trapped ions. The electronic configuration of indi- 
vidual ions, and the quantum state of their collective 
CM vibrations can be manipulated by coherent inter- 
actions of the ion with a laser beam, in a standing wave 
configuration, which can be pointed at any of the ions. 
The CM mode of axial vibrations may then be used 
as a "bus" to implement the quantum logical gates. 
Once the quantmn computation has been completed, 
the readout is performed through the mechanism of 
quantum jumps. Several features of this scheme have 
been demonstrated experimentally, mostly using a sin- 
gle trapped ion [||J|]. 

The unavoidable interaction of a quantum computer 
with its environment places considerable limitations on 
the capabilities of such devices 0. In this letter we 
make a quantitative assessment of these limitations for 
a computer based on the Cirac-ZoUer cold-trapped-ion 
design, in order to determine the best physical imple- 
mentation and the optimization parameters for quan- 
tum algorithms [||. There are two fundamentally dif- 
ferent types of decoherence during a computation: the 
intrinsic limitation imposed by spontaneous decay of 
the metastable states |e) of the ions; and practical 
limitations such as the random phase fluctuations of 
the laser driving the computational transitions or the 
heating of the ions vibrational motion. One could, in 
principle, expect that as experimental techniques are 
refined, the effects of these practical limitations may 
be reduced until the intrinsic limit of computational 
capability due to spontaneous emission is attained. 

The influence of spontaneous emission on a quantum 
computation with trapped ions depends on: the nat- 
ural lifetime of the |e) qubit; the number of ions, L, 
being used; and the quantum states of those ions. The 
number of ions which are not in their ground states 
varies as the calculation progresses, with ancillary ions 
being introduced and removed from the computation. 
The progression of the ions' states can be characterized 
well by an effective number of ions, Le, which have a 
non-zero population in the excited state |e). In the case 
of Shor's factoring algorithm ||^ (using long multiplica- 
tion), a reasonable estimate is « 2L/3. Therefore, 
to estimate the effect of decoherence during the im- 
plementation of Shor's algorithm, we will consider the 
following process: a series of laser pulses of appropriate 
strength and duration ( tt/2 pulses) is apphed to 2L/3 
ions, causing each of them to be excited into an equal 
superposition state {\g) + |e))/2. After an interval T, a 
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second series of laser pulses pulses) is applied, 

which, had there been no spontaneous emission, would 
cause each ion to be returned to its ground state. This 
is the "correct" result of our pseudo-computation. If 
there were spontaneous emission from one or more of 
the ions, then the ions would finish in some other, "in- 
correct" state. This process involves the sort of super- 
position states that will occur during a typical quan- 
tum computation, and so the analysis of decoherence 
effects in this procedure will give some insight into how 
such effects influence a real computation. A simple 
calculation shows that the probability of obtaining a 
correct result is 



P{T) w 1 - LT/Gto 



(1) 



where is the natural lifetime of the excited state |e). 
Thus the effective coherence time of the computer is 

The total time taken to complete a calculation will 
be approximately equal to the number of laser pulses 
required multiplied by the duration of each pulse. The 
time taken to switch the laser beam from ion to ion 
is assumed to be negligible. There are two types of 
laser pulse that are required in order to realize Cirac 
and ZoUer's scheme. The first requires pulses that are 
tuned precisely to the resonance frequency of the |e) to 
\g) transition, configured so that the ion lies at the node 
of the laser standing wave ( " V-pulses" ) ; the second re- 
quires pulses tuned to the CM phonon sideband of the 
transition, arranged so that the ion lies at the antinode 
of the standing wave ( "U-pulses" ) . The interaction 
of U-pulses with the ions is considerably weaker than 
the V-pulses, and so, assuming constant laser intensity, 
the U-pulse duration must be longer. Hence, in calcu- 
lating the total time required to perform a quantum 
computation, we will neglect the time required for the 
V-pulses. Because the entire calculation must be per- 
formed in a time less than the coherence time of the 
computer, we obtain the following inequality: 



(2) 



where is the total number of U-pulses, each of which 
has duration t^. The Hamiltonian for the interaction 
of these pulses with the ions is given by the following 
expression (ref. §, eq. (0)): 



H = 



hr] 
2y/L 



n[\e){g\ae~^'^ + \g){e\a^e^^] 



(3) 



In this formula, is the Rabi frequency for the laser- 
ion interaction, L is the number of ions in the trap, a 
(a^) is the annihilation (creation) operator for phonons 
of the CM mode and rj — ^JUlo^ cos^ d l2M(?v^ is the 
Lamb-Dicke parameter (here lo is the laser angular fre- 
quency, Q the angle between the laser and the trap 



axis, is the angular frequency of the ions' axial CM 
mode and M the mass of each ion). A careful calcu- 
lation, based on a perturbative analysis of the excita- 
tion of phonon modes other than the CM mode, shows 
that this Hamiltonian is valid if {Q,ri/2vx^fL)'^ <^ 1. 
The longest duration laser pulse that will be required 
to implement a quantum computing algorithm using a 
Cirac-Zoller quantum computer is a U-pulse of dura- 
tion tu = 2-1:^/1, l^lri. We will assume that all of the 
U-pulses required for the calculation are of this dura- 
tion. Therefore the lower bound on the duration of 
laser pulses is tu = yTr/vx, where y is a dimensionless 
"safety factor". This result can also be obtained from 
the naive uncertainty principle argument that there 
must not be appreciable power at the frequencies of 
the adjacent lattice vibrations. 

In order to attain the highest possible computational 
capability, one will need to minimize the duration of 
each laser pulse. Hence, it will be advantageous to 
employ an ion trap with the largest possible value of 
the trap frequency v^. However, the axial frequency 
cannot be made arbitrarily large because, in order to 
avoid crosstalk between adjacent ions, the minimum 
inter-ion spacing must be much larger than the size of 
the focal spot of the laser beam. The minimum separa- 
tion distance between two ions occurs at the center of 
the string of ions, which can be calculated by solving 
for the equilibrium positions of the ions numerically, 
resulting in the following expression : 
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where Z is the degree of ionization of the ions, e is the 
electron charge and Co is the permittivity of a vacuum. 
The spatial distribution of light in focal regions is well 
known [0. The approximate diameter of the focal 
spot is a^spot ~ Ai^, where A is the laser wavelength and 
F the f-number of the focusing system (i.e. the ratio 
of the focal length to the diameter of the exit pupil). 
Hence the requirement that the ion separation must be 
large enough to avoid cross-talk between ions, i.e. that 
3^inin 3> Xgpot, Icads to thc foUowing expression for the 
duration of the U-pulses: 



Z^ 



(5) 



where A is the atomic mass number of the ions. From 
eqs. I and I we obtain the following constraint on the 
number of ions L and the total number of U-pulses : 



7V„ii-84<2.0[s-im3/2] 
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We will now apply this bound to Shor's factor finding 
algorithm B . Let be the number of bits of the integer 
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we wish to factor. A careful analysis of the imple- 
mentation of the algorithm (using long multiplication) 
reveals that the required number of ions and U-pulses 
are given by: 

L = 51 + 2 , (7) 

Nu = 544^3 + 78P + lOl . (8) 

Note that there are asymptotically much more efhcient 
implementations, but they do not become competitive 
for the small number of binary digits under considera- 
tion here. If the measured Fourier transform and 
interleaving measurements were to be used in the com- 
putation, the number of ions required can be reduced 
to 4. However the intermediate measurements may 
increase the decoherence of the other ions due to scat- 
tered photons or unintended heating of the ions. It is 
for this reason that we have avoided use of this tech- 
nique in the assumptions underlying the algorithm. 

Equations |7| and |^ define a curve in (L, Nij) space, 
which taken in conjunction with the inequality (H) al- 
low us to determine the largest number of ions that can 
be used to implement Shor's algorithm in an ion trap 
computer with bounded loss of coherence. The linear 
relationship between L and eq.(|^), can then be used 
to determine the largest number that can be factored. 

As specific examples, we will consider the intrinsic 
computational capacity of Cirac-ZoUer quantum com- 
puters based on the following three ions: 

(i) Hg II: Z=l, A = 198; |e) is a sublevel of the 
Sd^Gs^ 2D5/2 state, \g) is the bd^^Gs'^ the two 
states being connected by an electric quadrupole tran- 
sition: A = 281.5 nm; Tq ~ 0.1s. 

(ii) Call: Z=l, A = 40; |e) is a sublevel of the 3d 2D5/2 
state, \g) is the 4s the two states being con- 
nected by an electric quadrupole transition: A = 729 
nm; t « 1.14s. 

(iii) Ba II: Z=l, A = 137; |e) is a sublevel of the 
5d ^D5/2 state, \g) is the 6s ^Si/2, the two states be- 
ing connected by an electric quadrupole transition: A 
= 1.76 /J,m; Tq ~ 47s. 

We shall assume that we have a very high numerical 
aperture focusing system, so that (although in practice 
such a high focal ratio would be difficult to achieve), 
and we will err on the side of optimism by putting 
the safety factor y — I. In figure 1 we have plotted 
the curves which limit the allowed values of L and , 
as given by eq.(^). We have also plotted, with a solid 
line, the "curve of factorization" defined by eqs. (^ 
and (||). The interception of the limiting curves for the 
different ions with the curve of factorization gives us 
the largest allowed value for the number of ions. Exam- 
ining these curves, we find that the size of the largest 
integer that can be factored by a Cirac-Zoller quantum 
computer based on Hg II, Ca II or Ba II ions is 6 bits, 9 



bits and 13 bits respectively. If these calculations were 
repeated with the less optimistic value for the safety 
factor, y = 3, one obtains 3 bits, 5 bits and 7 bits for 
the three species of ions, respectively. We note that 
the spontaneous emission lifetime is proportional to an 
odd power of A ( A^ in the case of electric quadrupole 
transitions) and so, for greater capability, eq.(^) sug- 
gests either going to longer wavelength (as seen with 
the three ion species above) or more highly forbidden 
transitions. 
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FIG. 1. The bounds on the numbers of ions, L, and the 
number of U-pulses, Nu , that may be used in a quantum 
computation without loss of coherence. The allowed values 
of Nu and L lie to the left of the curves. Curves for three 
ions are plotted. The unbroken line is the "factorization 
curve", specified by eqs. (7) and (8), which represents those 
values of L and which are required for execution of Shor's 
algorithm; the heavy black dots on this line represent the 
values of L and required to factor a number of 1 bits (1 — 
1, 2, ...15). 

As a more dramatic illustration of the theoreti- 
cal possibilities of the Cirac-Zoller scheme, one may 
consider a computer based on the 4f^^6s ^Si/2 ^ 
4f^^6s^ ^Fi/2 electric octupole transition of Yb II. This 
very long lived transition, which has received consid- 
erable attention because of its potential applications 
as an optical frequency standard, has a wavelength of 
467 nm and a calculated hfetime of 1533 days ]l^ . 
Performing a similar calculation to that given above 
suggests that, using this ion, it might be possible to 
factor a 438-bit number. Because such a calculation 
would require around 2200 trapped ions and 4.5 x 10^" 
U-pulses, taking about 100 hours, it would be difficult 
to over-emphasize the problems attendant on such an 
experiment. 

One may calculate the limits on factoring due to 
other causes of decoherence by a similar procedure to 
that used above. In this case, we will assume that the 
loss of quantum coherence due to sundry effects such as 
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random fluctuations of the laser phase or the heating 
of the ions' vibrational motion can be characterized by 
a single coherence time Tg. The effects of other causes 
of error, such as imprecise measurement of the areas of 
TT-pulses, which do not result in decoherence but never- 
theless lead to incorrect results in a computation, can 
also be characterized by the time Tg. Thus, in place 
of eq.(||) we now have the inequality . Using eq.(||) 
we obtain the following constraint on the values of the 
number of ions L and the number of laser pulses which 
can be used in a factoring experiment without signifi- 
cant loss of quantum coherence: 



0-84 <0.34[s-im3/2] 



y5/2^1/2^3/2;^3/2 
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Using the "factorization curve" specified by eqs. (7) 
and (8), one can obtain as before a value for the number 
of bits I in the largest number which may be factored. 
In this case the value of I will depend on the value of 
the coherence time Tg. In figure 2 we have plotted the 
values of / as a function of the experimental coherence 
time for the three species of ions discussed above. As Te 
increases, the largest number that can be factored also 
increases, until the limit due to spontaneous emission 
discussed above is attained. The slowest heating rate 
for a single trapped ion so far reported is 6 phonons 
per second (i.e. Tg = 0.17s) ||l^], and the laser phase 
coherence times longer than 10~^s have been achieved 
by several groups Q . Comparing these numbers with 
fig. 2, we see that, in principle, current technology is 
capable of producing a quantum computer that could 
factor at least small numbers (several bits). Note that, 
in contrast to the spontaneous emission bounds from 
eq. (6) (where To is, for quadrupole transitions, propor- 
tional to A^), eq.(9) argues for using shorter wavelength 
transitions. So we see from figs. 1 and 2 that Ca II is 
a good choice of ion for the experimental study of this 
technology because it allows a large number of opera- 
tions to be performed with realistic laser stability and 
ion heating requirements. 

The various causes of experimental decoherence 
which are mentioned above are all the subject of on- 
going research. It is not clear, for example, how laser 
phase fluctuations will affect quantum computations; 
it may be the case that the laser need be coherent 
only over the period required to execute each quantum 
gate operation. Furthermore, the heating rate of the 
ions' vibrational motion as a function of the number of 
trapped ions is not known. Other methods of coher- 
ent population transfer, which may be less susceptible 
to the effects of phase fluctuations, for example stim- 
ulated Raman adiabatic passage (STIRAP) jl^, are 
being investigated. 
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FIG. 2. Tlie variation of the number of bits 1 in the 
largest integer that may be factored with the experimen- 
tal coherence time for the three ions discussed in the text. 
The maximum values of the computational capacities for 
the ions Hg II and Ca II are the limits determined by spon- 
taneous emission. 

It is clear that if quantum computation is to over- 
come decoherence and other errors, then error correc- 
tion must be used extensively. So far, suggestions for 
error correction have relied either on variations of the 
"watch dog" effect or on exploiting the prop- 

erties of certain entangled states to reduce the impact 
of decoherence in a quantum memory ||l^. The lat- 
ter has not yet proven to be practical for use during 
a quantum computation, primarily because there has 
not been any analysis of the success of the method 
under realistic assumptions on operator errors. If op- 
erational errors were negligible, the effect of decoher- 
ence on quantum memories could be reduced arbitrar- 
ily. However some of the "watch dog" methods that 
have been suggested are quite practical. For example, 
many computations require the use of ancillary qubits 
which are periodically returned to the ground state. 
Measuring these ancillas when they are supposed to be 
in the ground state can be used to help dissipate er- 
rors. Recent simulations indicate that this method 
is indeed helpful in maintaining the state of the com- 
putation. Implementation of the method does require 
intermediate measurements. In any case the effect of 
using a "watch dog" method is to stabilize the effec- 
tive decoherence time Tg by ensuring less dependence 
between the errors of successive operations. 

In conclusion, we have derived quantitative bounds 
which show how the computational capabilities of a 
trapped ion quantum computer depend on the relevant 
physical parameters and determine the computational 
"space" [L) and "time" {Njj) combination that should 
be optimized for the most effective algorithms. The 
effect of this bound has been illustrated by calculating 
the size of the largest number that may be factored us- 
ing a computer based on various species of ion. Our re- 
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suits show there is reason for cautious optimism about 
the possibihty of factoring at least small numbers using 
a first generation quantum computer design based on 
cold trapped ions. However, the large number of pre- 
cise laser operations required and the number of ions 
involved indicates that even this computationally mod- 
est goal will be extremely challenging experimentally. 
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